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After the previous work on gravitational frequency shift, hght deflection (arXiv: 1003.5296) and 
perihelion advance (arXiv:0812.2332), we calculate carefully the fourth gravity test, i.e. radar echo 
delay in a central gravity field surrounded by static free quintessence matter, in this paper. Through 
the Lagrangian method, we find the influence of the quintessence matter on the time delay of null 
^— ^ , particle is presence by means of an additional integral term. When the quintessence fleld vanishes, 

vN I it reduces to the usual Schwarzschild case naturally. Meanwhile, we also use the data of the Viking 
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lander from the Mars and Cassini spacecraft to Saturn to constrain the quintessence field. For the 
Viking case, the field parameter a is under the order of 10"". However, a is under 10~^* for the 
Cassini case. 
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I. INTRODUCTION 



*Zi , The discovery of dark energy is a landmark of cosmology. In 1998 jl|, it became known that our universe is ac- 

celerating according to the observation data from the distant supernova. This fact was proved later by further la 
(•^ ■ Supernovae (SNe la) y and microwave background (CMB) [3| of WMAP. The special matter contained a negative 
l/~J ' pressure which can drive the universe acceleration, called dark energy. The large scale distribution of galaxies obser- 
f^ ' vations [4| also tells us there is about 70% dark energy in the universe, which seems more like an evenly distributed 
background component and it cannot assemble in the large scale. Else, the theories of big bang nucleosynthesis and 
galaxy formation confine in the early universe the ratio of dark energy to very small values and only it becomes 
large after galaxy formation. In the theoretical analysis, the simplest candidate is the cosmological constant. If so, 
our universe will be further accelerating and the whole space will more like a de Sitter topological structure. The 
cosmological constant combined to cold dark matter leads to a so-called ACDM cosmological model 5] . In spite of the 
success of ACDM model in many respects, there are still two big troubles in this model; i.e. the famous fine-tuning 
and cosmic coincidence problems. Hence, many theoretical physicists prefer a dynamical dark energy model to a 
ACDM model. The simplest dynamical model is a quintessence field which is a scalar field with canonical momentum 
y. When quintessence rolls along itself the potential curve, it will provide a negative pressure if the potential curve is 
very flat i.e. the quintessence field is at a slow roll state. Else, its attractor solution can solve the cosmic coincidence 
very well. 
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As we know that, the cosmological constant model can give a static sphericahy symmetric black hole i.e. 
Schwarzschild-de Sitter black hole (SdS-BH) which is corresponding to a usual Schwarzschild case. However, how 
to obtain a counterpoint to that of Schwarzschild case just like the SdS-BH in cosmological constant model? After 
this problem occurred some fundamental work had been done by people such as Gonzalez 7|, Chernin et.al.[8| and 
Kiselev [9| and so on. In the early works, there were no horizon and no 'hair' problems [7|, lg| which means no black 
hole could be embraced into these metrics. Later, in order to solve this problem Kiselev employed nonzero off-diagonal 
energy momentum tensor in which its coefficients satisfied additivity and linearity conditions. Hence, a steady black 
hole solution is obtained [9|]. The process is introduced simply in part H. 

Since Kiselev [9[ obtained the exact black hole solutions encoded in quintessence, many people have performed 



studies on its special properties including: quasinormal frequencies fONMl llOl, th ermodynamics entropy [111 ], gravity 



geodesic precession [12| and the classical three gravity tests on solar system [13|,|14| and so on. The last phenomenology 
aspects refer to the gravitational frequency shift, the deflection of light and the precession of perihelion of Mercury. 
In Ref.[l3| via calculating the photon's 4D momentum, it is found that the photon frequency increases for the range 
LUq e [—1,-1/3] and decreases for the range [—1/3,0], where tUq is the quintessence state parameter satisfying the 
equation of sate (EOS). The field parameter a is sensitively dependent on the state parameter ujq. Comparing with 
the H-masers of the GP-A redshift experiment data, the constraint on quintessence matter is a G (10^^^, 10^^) with 
the range ujq G [—1,0]. Meanwhile, for the light deflect test its influence behaviors also are different according to 
various ujq values. The analytic results are obtained via the integral photon's Binet equation with special solvable 
value ujq = {—1, —2/3, —1/3, 0}. By using the long-baseline radio interferometry data, the corresponding constraints 
on the quintessence field parameter a = { — , lO^'^^, 10^"^, 10^^} are obtained where the label "— " means that there 



is no infiuence on the defiection of light by quintessence for the case ujq = ~1 13|. Then after that, the test of 



perihelion advance has been presented in Ref. [14]. The trajectory of the test mass and the motion of a binary system 
are discussed in detail. However, as far as we know there is no work relating to "the fourth test" [ISl . Ilq in gravity 
field i.e. time delay in gravity test. Considering the above various factors, we calculate carefully the behavior of 
gravity time delay infiuenced by quintessence matter and try to find the constraints on quintessence field parameter 
in alternative view here. 

This paper is organized as follows: in section II, we present the Kiselev solution briefly. In section HI, we use the 
Lagrangian method to find the analytic expression of time delay in different points. In section IV, we use the data 
of the Viking lander on Mars and the Cassini spacecraft on the way to Saturn to constrain the quintessence field 
parameter. Section V is the conclusion. We adopt the signature (-I-, — , — , — ) and put h, c, and G equal to unity. 

II. THE STATIC SPHERICALLY SYMMETRIC BLACK HOLE SOLUTIONS WITH FREE 

QUINTESSENCE MATTER 



The modified Schwarzschild space surrounded by static spherically symmetric quintessence matter is given by 
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where M is the black hole mass, LOq is the state parameter of quintessence, a is the corresponding quintessence field 
parameter. In this paper, we evaluate the quintessence field with cjq in the range of [—1,0] without considering the 



supernova dimming. Interestingly, this metric reduces to the pure Schwarzschild space and the Schwarzschild-de Sitter 
space for the cases of a = and ujq = —1, respectively. The energy momentum coupling with evenly distributed 
quintessence is shown as 



Tl = P,{r), 
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where pq is the density of quintessence matter. The internal structure B function satisfied B = — (Sojg 
the additivity and linearity principle T/ = T^ . The isotropic average over the angle components is 



<T > 
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where the relationship < Vir^ >— ■nSfrnr''^ is used and the state equation is Pq — ujqPq where ujq — 7/8. 
With the relationship X + v — 0, the energy momentum tensor components are 



r/ = r; = p,; t,« = r; = --p,(3^, + i). 



The quintessence density is 



a 3w„ 



Pq = o: 
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So if we require the density of energy positive, pq > 0, we deduce that a is negative for u)q negative. The curvature 
has the form of 

1 — 'Sujn 






(7) 



Apparently, r = is the singularity for a;^ 7^ { — 1, 0, 1/3}. The free quintessence creates an outer horizon of de 



Sitter universe at r = r^ for 



-l<OJq<--, 

and also generates an inner horizon of black hole at r — Vq for 

- - < Wg < 0. 



(8) 



(9) 



III. RADAR ECHO DELAY ON SOLAR SYSTEM 



Except for the gravitational redshift 



Q 



13| , perihelion procession [lj| and bending of light [13| , a further measurable 



effect concerning the null geodesies is the "time delay" of radar signals. Unlike the flat space, the travel time of light 
between any two given points increases due to the presence of space curvature. This type of travel increasing time 
can be measured via using the propagation of radar signals in solar system, which is proposed firstly by Shapiro in 
1964 !l5|. Else, it is also constituted by the "fourth test" of GR known as radar ranging [16|. Here, the time delay 
for radar ranging is calculated in the space of metric ^ by the Lagrangian method. The Lagrangian for a test mass 
particle in the field described by the equation ([T]) is 

1/2 
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where the overdot represents differentiation with respect to the afBne parameter ^ along the geodesies. The Euler- 
Lagrange equations read as foUows 

±(dCA_dC^^Q (11) 

Without loss of generality, the observer is confined to orbits with 6 ~ 7r/2 and = 0. Hence, the Lagrangian Cp 
becomes 

Cp=iy-f{r)i+^^+r'^^ . (12) 

According to the cyclic coordinates ^p and t, we can identify two constants of motion, 

r'^^L. (13, 

2M a ,dt 

The third motion equation can be given by the normalization relation of photons g^uP^p'^ — 0, 

drV _ 2 L"^ ( 2A/ a \ 

dc) -^ '^[^' — 'v^^;+i)- ^^^> 

The dynamic evolution of photon is determined fully by Egs. p^ . ([T?l) and P^ in this space. Combining Eq. p^ 
with (jlSp . the photon's trajectory equation can be gotten as 

1 dry _ fEV 1 / 2Af a \ 

Meanwhile, we define two new parameters, one is the impact parameter b = L/E, which means the effective sighting 
range, the other is the photon's effective potential 1/ B'^ir) where 

^« = '^(i- — -;3^3;ttJ • (17) 

According to Eq. (|17l) . the photons' orbital equation p^ can be rewritten as 



The distance of closest approach rg is defined by 

So we can get the relationship h — B{ro). According to Eqs.(fT4|) and (fT5|). the differential relation between time t and 
radius r is obtained as follows: 

where /(r) = 1 — 2M/r — Q;/r^"''+^. So the time travel for the null geodesic is 

dr 
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The first term t{r,rQ)GR is the usual case in GR 18] and the last one is the effect of quintessence field which will 
be solved particular in the remained parts. Apparently, once the quintessence field disappears, the result is reduced 
to the Schwarzschild one. Because the undetermined state parameter ujq which is coupling in the exponent of r or 
To, we select the solvable integers 3ujg + 1 = (—2, —1, 0, 1) with Wg € [~l7 0] to integrating Eq. (|22p . This kind of 
method, which also is shown in Ref. 13| and 14|, is useful to solve approximately sophisticated differential equation 
in which the independent variable contains some unknown parameter. The integrating results are shown in TABLE ID 
The corresponding sketchy outline of the additional term i(r, ro)QM versus radial r is shown in Fig. ([l}. It illustrates 
clearly that the influence of quintessence on time delay heightens with the decreasing quintessence state parameter 
LUq correspondingly. 

TABLE I: The forms of t{r,ro)QM 
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IV. CONSTRAINTS FROM EXPERIMENT OF VIKING LANDER ON MARS 



The time interval between emission and return receiver is 



AT = 2i(re,ro) + 2i(rfl,ro) 
- ATGii + Ar„, 



(23) 
(24) 



where the first term ATga is the usual GR value the modified term ATq depends on the state parameter LOq of 
quintessence which is presented in TABLE HIl with the limits of r® 3> tq and r^j 3> tq. The proper time Ar which 



elapses on Earth is in relation to the change in coordinate time AT in the form Ar = v/5oo(^©)Ar. 

So the excess round trip Ar for signal emitted from Earth and bounced off the other planet can be obtained as 



Ar = ATGfl(l + ATD) 



(25) 
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FIG. 1: Additional term t{r,ro)QM versus radial r coordinates for a = 10 '^ and solar radius ro — 1. 



TABLE II: The estimates on quintessence field parameter a from Mars 
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|a| < 3.78879 x lO-^^ 



The GR value Atgr is 
Atgr = 4M 
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The best experimental constraints from Solar system on the time delay is the Viking lander on Mars and give 
\Atd\ < 0.002 [17]. The distances of Earth and Mars from Sun are r© = 1.525 x lO^^cm and tr = 2.491 x IO^^cto, 
and the Sun radius is ro — 6.955 x lO^'^cm. The final additional time delay AT^ caused by quintessence field is shown 
in Table ini Meanwhile, we also give the corresponding constraints on field parameter a in TABLE |ll] based on the 
data of Viking lander on Mars [iT^. Else, we can read that the smaller the parameter ujq, the stronger constraints on 



V. CONSTRAINTS FROM CASSINI SPACECRAFT EXPERIMENT 

It is also known that the highest precision data actually have come from the tracking of spacecraft which emit 
signal rather than reflection of radar off of planets for its surface influence [171, |l8| . Here we use the data of Cassini 
spacecraft from Cassini-Huygens mission which are composed of two main elements: the NASA Cassini orbiter and 
the ESA Huygens probe. The changed frequencies for Cassini spacecraft experiment are [19| 



" v[Q) dt 
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where :/(0) is the emitted frequency from Earth and i'(i) is the received frequency by Earth, respectively. The 
contribution to frequencies by quintessence, i.e. a-term, is 

d 



Vo 



dt 



-AT„ 



(28) 



It is known that the time delay experiment is measured nearly 25 days i.e. 12 days before conjunction and 12 days 
after. The Cassini, the Sun and the Earth are almost aligned at conjunction and the geocentric distance of spacecraft 
is 8.43 AU, with a minimum impact parameter bmin = l.Bi?©. According to these measures, the general relativity 
(A = 1) contribution to y{i') is 6 x 10~^° within an accuracy of 10"^"*. This leads immediately to a upper bound 



2/a(12rf)-y„(0)<10- 



(29) 



Because of the spacecraft being much farther away from Sun than Earth, we adopt that t'q — dr^/dt w uq w 
29.78 Km/ sec is approximately Earth's orbiting velocity around Sun. For the distance of Cassini is adopted approx- 
imately its conjunction position i.e. rji — %A2>AU — \AU = 1.1 x lO^^m. According to the various quintessence state 
parameters ujq and the above various conditions, four types of additional frequency y^ are listed in Table III. Except 
for the case of ujq — 0, the trends of constraints on a versus ujq are similar with the experiment of the Viking lander 
on Mars. 



TABLE III: The estimates on quintessence field parameter a from Cassini spacecraft experiment 
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VI. CONCLUSION 



In this paper, we have used the Lagrangian formulation to analyze the influence of quintessence field on "the fourth 
gravity test" i.e. radar echo delay in a central gravity field without considering the internal charge. By using the data 
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of Viking lander and Cassini mission, we have obtained the constraints of quintessence field from alternative channels. 
Now, two points need to be emphasized especially as follows. 

1. There are two key parameters Uq and a in this model. The formal Uq, which connects the pressure pq and 
energy density pq via the state equation pq — itJqPq, is the state parameter of the quintessence. The latter a is the 
quintessence field parameter, which indicates the influence of quintessence matter over the space outside the black 
hole. Here we need to pay more attention on the cause and effect of a. In Kiselev's original work [9|, one additively 
and linearity condition T/ = T^ is adopted to solve the previous "no hair no horizon" problem [71, |8J . Combining the 
Einstein equation in spherically symmetric static space, the quintessence energy momentum tensors 1^ and ([3]) are 
finally reduced into a differential equation. 






^'3:^ + 3(l+^g)'^^ + (3^g + l)/ = 0, (30) 

where function / is decided by metric function e^ via A = — ln(l + /). Then Eq. (j30p has two solutions, 
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fi ^ ^3w„+i ' v^-*-) 
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Here, a and Vg are the normalization factors. Vg is the usual Schwarzschild solution and a represents the quintessence 
field parameter. The expression (j3ip is the original position where parameter a appear firstly in this model. It is helpful 
to understand its meaning from two aspects. One is the density of the quintessence energy pq — 3Q!ajg/2r^(^+"''-'. So 
the sign of normalization constant a and the matter state parameter Uq have to satisfy the condition of aujq ^ 0. 
Considering the fixed state parameter, the quintessence increases with bigger a. Hence, its magnitude manifests 
enough the strength of quintessence field. The other feature is the curvature of space R = 3aa;g(l — 3a;g)/r^^"''+^^ If 
one fixes the state parameter, we can find the curvature R will be larger with bigger \a\ where this point also means 
the space will be bended more. From the expressions of density we can find that the dimension of a is sensitively 
dependent on the value of state parameter ujq. The specific dimensional form of a is expressed in Table 6.1 in Ref. 
[13j according to various state parameters. If the quintessence field could be treated as the dark energy and the 
critical density of the universe could be assumed as the value of pq, we could find that with bigger astronomical scales 
the magnitude of |a| is larger. In another words, the quintessence field can have infiuence more deeply on the large 
astronomical scale. 

2. In this model, we only consider a neutral black hole as well as SdS-BH. The internal charge is not involved in this 
paper. If we want to consider the charged case, it should refer to the general solution for the Reissner-Nordstrom-de 
Sitter black hole surrounded by the quintessence matter [9|, which is given by 

r r^ a^ \ r . 

where a? = 3/A. If we choose the simplest singular Schwarzschild case [13|, |lj|, the last term can be reduced to 
(^) ' — > 3^q+i just like this model. This solution ([55]) also has a rich spectrum of limits. If e — > 0, there will 
be no charge. If a? — >■ oo, there will be no de Sitter curvature. If Vg — > 0, e — >■ 0, it will become a self-gravitating 
quintessence without black hole. Hence, if we consider the internal charge, we must start from another metric space, 

g^^^l-^ + -~-^-^. (34) 



.r^l-^ + i-^-f^)- , (33) 
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